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Angular asymmetries are simple, intuitive, model-independent observ- 
ables used to identify spins of new elementary particles. In the case of 
Drell-Yan-like boson resonances, we generalize the well-known center-edge 
4^ ' angular asymmetry to optimize spin identification when only a limited sam- 

Q-i pie of events is available. By choosing simple weight functions W{&) in 

^ • integrals over the polar angle 8, such as W — cos™ 9, we can improve spin 

_ discrimination significantly in production and decays of spin 0, 1, and 2 

. bosons. The power n can be tuned in particular cases, but n = 2 (n = 1) 

works well for any forward-backward symmetric (non-symmetric) decay to 
massless particles. 

New particles, such as Higgs, Z', excited gravitons, sparticles, etc., may 
well be discovered in collider data from pp (Tevatron), pp (LHC), or e + e~ 
(ILC) scattering. The spin of a new virtual particle is often determined in 
the experiment from the analysis of angular distributions in the heavy parti- 
cle's decay. This analysis follows a standard quantum-mechanical expansion 
over spherical harmonics that depends on the spin of the massive particle 
state. In this contribution, we describe an analysis procedure for decay 
angular distributions that is more efficient in using limited event statistics 
than the commonly adopted method. One way to determine the spin is to 
exploit the "center-edge" asymmetry that typically exists between forward 
and transverse decay products from a boson [TJ [2] . We show how to modify 
the usual "central-edge" asymmetry definition to increase its discriminating 
power for the spin. We keep the discussion transparent by working at the 
parton level, but consider a simple analytic model for typical experimental 
acceptance. We checked that our main conclusions remain valid in a more 
realistic calculation, by comparing our method to a fully differential calcu- 
lation using ResBos [3] that includes proton PDFs, NLO corrections, and 
NNLL resummations in QCD. 
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For decay of a boson into two back-to-back massless particles, we define 
a general asymmetry by an integral 

A = ^ W(z)P(z)dz. (1) 

Here z = cos 6, 9 is the decay polar angle in the boson center-of-mass frame, 
P(z) = is the normalized production probability density, and W(z) is 
a weight function. One popular choice for W(z) has been the center-edge 
step function, [H m, H El El El E] 

W CE (z)= +1 (|*| >z*); -1 (|z|<z*), (2) 

with z* ~ 0.5. We will propose an alternative definition for W(z) that is 
more optimal than Wce( z )- 

To discriminate between two possible boson spins a and b in a measure- 
ment, we introduce respective probability densities P a (z) and Pb(z) and 
their asymmetries A a and A&, and construct the ratio 

Rab = SA a+ ti b (3) 
The expected statistical uncertainty is estimated by 



SA = J J P(z) W(z) 2 dz - ^ P(z) W(z) dz^j . (4) 

R a b is independent of additive and multiplicative constants in W{z). A 
larger value indicates that A with the chosen weight is more sensitive 
to spin. We wish to choose W{z) to maximize R a b, but we don't want 
to fine-tune W(z) for each particular process/model/experiment. We will 
show for typical cross-sections of interest that the simple choice = z n 

works well. Although n may be tuned, n = 2 (n = 1) is largely sufficient for 
distributions that are forward-backward (non) symmetric under z — > — z. 

A choice of W(z) that amplifies in the range of z where the difference 
AP(z) = Pb — Pa is large in magnitude can increase the numerator in eq. Q. 
(Note that AP is normalized to zero, so it must change sign). However, 
amplifying this range of z will unbalance the cancellation of statistical fluc- 
tuations across the whole —1 < z < 1 range and increase the denominator 
in eq. ([3]), which must happen since AP is normalized. On the other hand, 
a uniform choice, such as Wqe in eq.©, gives equal weight to the whole 
z range, reducing the denominator through cancellation of fluctuations but 
not optimimizing the numerator. It appears there is a minimization problem 
to solve. 
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Let us consider a concrete 
example. We compare the low- 
est order QCD densities Pq 
and Pi for spin-0 and spin- 
1 bosons produced by massless 
partons and decaying to mass- 
less fermion pairs: 

p„ 1 



Pi 



2' 

|(l + ^ 2 + 2 C 




d 



c + d 



z). 



The term arises for general 
boson-fermion chiral couplings, 
e.g., it is nonzero in a parity- 
violating decay. Consider first 
the symmetric case c = d, so 
that Pi = (3/8) • (1 + z 2 ). We 
try a class of weight functions 



Fig. 1. The solid line is the statistical signifi- 
cance ratio i?oi for the particular case of spin 
and 1 bosons decaying to forward-backward 
symmetric fermions. i?oi is plotted as a func- 
tion of n in the weight functions W = \z\ n . 
The points are obtained for discrete choices of 
n after modulating the leading-order densities 
Pa,b by a function ((?]) simulating experimental 
acceptance, for a = 6 and /? = 0.2. 



W{z) 



(5) 



with n selected to emphasize contributions from the z intervals with large 
AP(z). Figure [JJ confirms the existence of an optimal weight that maximizes 
Poi- Recalling that confidence limits are determined by Rs/N in the case 
of N sample events, Wqe would require about 1/3 more events to achieve 
the same level of significance as the optimal weight - see Table [TJ A similar 
conclusion holds for parity violating processes (c ^ d), with the results for 
the maximal symmetry violation d = shown in Table [TJ [In this case, we 
modified the center-edge weights as W(|z|) — > sign(z) • W(|^|)]. Provided 
there is a significant monotonic variation of AP from the center to the edge 
of the range of z, the optimum weight will be close to quadratic for c = d 
and linear for d = 0. 

These conclusions stay valid when the above lowest-order numerical esti- 
mates for R a b are modified by higher-order QCD corrections to P(z), smear- 
ing by parton distributions in the case of hadronic collisions, and detector 
acceptance constraints. For example, the angular dependence of P(z) can 
be modified by the limits y m \ n < y < y max (z) on the boson rapidity. We can 
explore the impact of these corrections by multiplicatively modulating the 
lowest-order parton-level z-dependence of the differential cross section by a 
function 

rM+AM roc 

m{z) = dM dy 6(y max - y) X Q(y - y min )^(y, M) (6) 

JM-AM J-oo 
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B^ff 


R a b with Wc£ 


Optimal n, R a b(n) 


% more events with Wqe 


a 


b 











1 


0.095 


1.6, 0.1097 


34% 





2* 


0.106 


4, 0.155 


113% 


1 


2* 


0.205 


3, 0.264 


66% 







0.45 


0.8, 0.52 


33% 


B — > 77 











2* 


0.239 


1.3, 0.272 


30% 


*(e q = 0.1) and **(d = 0). 



Table 1. The R a b ratio for a boson B decaying to difcrmions // or diphotons 77, 
computed for various spins a and b, using Wce and W = \z\ n weights with the 
optimal n value. Also shown is the percentage increase in the number of events 
that would be needed when using Wce to achieve the same statistical significance 
as with the optimal W — \z\ n . 



with a boson mass M, AM bin size around the resonance peak, step func- 
tions Q(y) indicating acceptance constraints imposed on the rapidity, and 
a function K(y, M) containing all details of QCD corrections, parton dis- 
tributions, boson propagator, etc. For a typical z distribution, m(z) can be 
modelled by 

m(z)oc(l-z a f, (7) 

where a and (3 depend on details of the calculation. Fig. [2] compares rep- 
resentative angular distributions for fermionic decays of massive bosons ob- 
tained with ResBos for a particular parameter set, which realistically in- 
cludes the effect of PDFs, QCD NLO and NLL corrections, and acceptance 
cuts, against calculations using the modulation model ([7]). We see that the 
m(z) function with a and f3 fitted to the ResBos points captures the general 
features of the ResBos predictions. For our purposes of comparing the statis- 
tical performance of different polar weights, the precise values of the powers 
in ([7]) are not found to be crucial. The effect of this modulation on the po- 
lar asymmetry is also shown in Fig. [TJ Although the statistical significance 
is reduced generally, as might be expected because the acceptance reduces 
the center-edge asymmetry, the conclusions about the optimum weight are 
essentially unchanged. 

To demonstrate limitations of the above conclusions in relation to the 
overall amount of center-edge symmetry, we may look at the decay distri- 
bution of a spin-2 boson to fermions, 

P2 = \(l-?>z 2 + Az i )e q + l(l- Z i )e g . (8) 
This depends in addition on the fractions of Drell-Yan events e q and e g 
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Fig. 2. Realistic P{z) for spin-0 (left graph) and spin-1 (right graphs) bosons 
decaying to massless fermions. The data points were generated using ResBos [3 
for boson mass M = 1 TcV produced in pp collisions at ^/s = 14 TeV. This 
includes NLO and NLL QCD corrections, CTEQ 6.6 PDFs [5], and &p T > 25 GeV 
cut on the fermion's transverse momentum. The solid lines represent a fit using a 
modulated P(z) with m(z) from eq. (J7]). 



produced via qq and gluon-gluon fusion respectively. Under the constraint 
e q + e g = 1 at a pp collider, e q typically varies from e q ~ at very low boson 
masses <C 1 TeV to e q ~ \ at 4 TeV. We find that the statistical significance 
of an optimized power law W = \z\ n is generally much better than Wqe 
across this range (see for example Table [T]). However, the advantage reduces 
as e q grows. This is because at large e q there is no longer any center-edge 
asymmetry at all and center-edge type observables are generally less useful. 
The optimal power is of course dependent on e q , but a simple universal 
quadratic choice W = z 2 works very nearly as well. 

Another important decay mode for identifying boson resonances is to two 
prompt photons, which are relatively clean to identify. The corresponding 
tree level probability densities in this case are 

= ~ Pi = 0, P 2 = ^(1 + 6z 2 + z')e g + - z*)e q . (9) 

The statistical significance of the optimized power law and Wqe choices for 
spin 2 vs. spin with acceptance modeled by (J7]) follows a pattern similar 
to the di-fermion final state (see Table [T|) . 

Aside from the simple powers, we have also investigated other possible 
angular weights, particularly for those that strongly enhance the regions 
where AP is large, and find they are not superior in statistical significance. 
Choosing a weight function that was orthogonal to one of the two raw prob- 
ability distributions used in any comparison did not produce better results 
for the R a b ratio when the experimental acceptance was taken into account; 
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part of the reason is that the latter destroys orthogonality. To conclude, the 
sensitivity of the z = cos distribution to the spin of the decaying heavy bo- 
son is increased by generalizing the center-edge asymmetry via eq. (H|) with 
the weight function W(z) = z n and n = 2 for a forward-backward sym- 
metric (n = 1 for a forward-backward non-symmetric) distribution. These 
findings are best noticed in a leading-order analysis of the decay angle de- 
pendence, but they are robust against radiative and acceptance corrections 
arising in the full NLO differential calculation. 
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